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Abstract 

The spontaneous generation of the chromomagnetic field at high temperature is investigated 
in a lattice formulation of the S'L r (2)-gluodynamics. The procedure of studying this phenomenon 
is developed. The Monte Carlo simulations of the free energy on the lattices 2 x 8 3 , 2 x 16 3 and 
4 x 8 3 at various temperatures are carried out. The creation of the field is indicated by means of 

IT) . 

, the x -analysis of the data set accumulating 5-10 millions MC configurations. A comparison with 

the results of other approaches is done. 

o 

O ' 1 Introduction 

Among interesting problems of modern cosmology the origin of large-scale magnetic fields is intensively 
attacked nowadays. Various mechanisms of the field generation at different stages of the universe 
■ evolution were proposed [TJ. Basically they are grounded on the idea of Fermi, Chandrasekhar and 

Zel'dovich that to have the present day galaxy magnetic fields of order ~ lfiG correlated on a scale 
^ ■ ~ IMpc seed magnetic fields must be present in the early universe. These fields had been frozen in a 

cosmic plasma and then amplified by some of the mechanisms of the field amplification. One of the 
ways to produce seed fields is a spontaneous vacuum magnetization at high temperature T [21 13 HI [5] . 
Actually, this is an extension of the Savvidy model for the QCD vacuum [6], proposed already at 
T = and describing the creation of the Abelian chromomagnetic fields due to a vacuum polarization, 
in case of nonzero temperature. At zero temperature this field configuration is unstable because of 
the tachyonic mode in the gluon spectrum. At T ^ 0, the possibility of having strong temperature- 
dependent and stable magnetic fields was discovered [4]. The field stabilization is ensured by the 
temperature and field dependent gluon magnetic mass. 

Another related field of interest is the deconfinement phase of QCD. As it was realized recently, 
this is not the gas of free quarks and gluons, but a complicate interacting system of them. This 
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was discovered at RHIC experiments [7j and observed in either perturbative [4 S 8] or nonperturbative 
[9] investigations of the vacuum state with magnetic fields at high temperature. In Refs. [4j [8] 
the spontaneous creation of the chromomagnetic fields of order gB ~ g 4 T 2 was observed in SU(2)- 
and SU (3)-gluodynamics within the one-loop plus daisy resummation accounted for. In Ref. [9] the 
chromomagnetic condensate of same order was obtained in stochastic QCD vacuum model and method 
of dimensional reduction by comparison with lattice data. In Refs. [10J the response of the vacuum 
to the influence of strong external fields at different temperatures has been investigated and it was 
shown that the confinement is restored by increasing the strength of the applied field. These results 
stimulated the present investigation. 

We are going to determine the spontaneous creation of magnetic fields in lattice simulations of 
SU (2)-gluodynamics. In contrast to the problems in the external field, in the case of interest the field 
strength is a dynamical variable which values at different temperatures have to be determined by means 
of the minimization of the free energy. This procedure is not a simple one as in continuum because 
the field strength on a lattice is quantized. To deal with this peculiarity, we consider magnetic fluxes 
on a lattice as the main objects to be investigated. The fluxes take continuous values, and therefore 
the minimization of the free energy in presence of magnetic field can be fulfilled in a usual way. These 
speculations serve as an explanation of the strategy of our calculations. 

One of the methods to introduce a magnetic flux on a lattice is to use the twisted boundary 
conditions (t.b.c.) [IT] . In this approach the flux is a continuous quantity. So, in what follows we 
consider the free energy F(tp) with the magnetic flux ipona lattice in the S'f7(2)-gluodynamics and 
calculate its values at different temperatures by means of Monte Carlo (MC) simulations. We will 
show that the global minimum of F(<p) is located at some non-zero value (p m in dependent on the 
temperature. It means the spontaneous creation of the temperature-dependent magnetic fields in the 
deconfinement phase. 

The paper is organized as follows. In sect. 2 some necessary information about the magnetic 
fluxes on a lattice is adduced. In sect. 3 the calculation details and the results are given. Section 4 is 
devoted to discussion. 

2 Magnetic fields on a lattice 

In perturbation theory, the value of the macroscopic (classical) magnetic field generated inside a system 
is determined by the minimization of the free energy functional. The interaction with the classical 
field is introduced by splitting the gauge field potential in two parts: A^ = A^ + A^}, where A 1 ^ 
describes a radiation field and A^ = (0, 0, Hx 1 , 0) corresponds to the constant magnetic field directed 
along the third axis. However, on a lattice, the direct detection of the spontaneously generated field 
strength by straightforward analysis of the configurations, which are produced in the MC simulations, 
seems to be problematic. Therefore, it is reasonable to follow the approach used in the continuum 
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field theory. 

First, let us write down the free energy density, 

Z{ V )= j[DU{ v )\^f{-S{V( v ))}. (2) 

Here, Z(ip) and Z(0) are the partition function at finite and zero magnetic fluxes, respectively; the 
link variable U is the lattice analogue of the potential A^. 

The free energy density relates to the effective action as follows, 

Fdp) = S{ip) - 5(0), (3) 

where S(<p) and 5(0) are the effective lattice actions with and without magnetic field, correspondingly. 

To detect the spontaneous creation of the field it is necessary to show that the free energy density 
has the global minimum at a non-zero magnetic flux, tp m i n ^ 0. 

In what follows, we use the hypercubic lattice L t x (L t < L s ) with the hypertorus geometry; 
L t and L s are the temporal and the spatial sizes of the lattice, respectively. In the limit of L s — » oo 
the temporal size L t is related to physical temperature. The one-plaquette action of the SU (2) lattice 
gauge theory can be written as 



l--TvU^(x) 



(4) 



U^ix) = U„(x)U u (x + ap)Ul(x + av)Ut{x), (5) 

where (3 = 4/c/ 2 is the lattice coupling constant, g is the bare coupling, U fi (x) is the link variable 
located on the link leaving the lattice site x in the fi direction, U^ v (x) is the ordered product of the 
link variables. 

The effective action 5 in |(3]) is the Wilson action averaged over the Boltzmann configurations, 
produced in the MC simulations. 

The lattice variable U^{x) can be decomposed in terms of the unity, /, and Pauli, aj, matrices in 
the color space, 



n, s / U°(x) +iUl{x) Ultx) + iUlix) 

U il {x)=IUl{x)+ici j U^x)=\ »l J £ ' 

-U*(x)+iUfa) U°(x)-iU*(x) 



(6) 



The four components U^(x) are subjected to the normalization condition - U^(x)U^(x) = 1. Hence, 
only three components are independent. 

Since the spontaneously generated magnetic field is to be the Abelian one, the Abelian parametriza- 
tion of the lattice variables is used to introduce the magnetic field, 

r(r) i cos<^(x)e i8 " (:c) Bin<f> li (x)e ix »W ^ 
' - sin0 M (x)e~ lx ^^' cos 4> tl (x)e^ l9 f 1 ^ 
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U„({x ,x„L„x,}+v) = U,({x ,x„O,x,}+v)e 

hm — identified ►<•- 



V 



U v (x l ,,x 1 ,L„x s ) 



identified 

) = u„(x ,x„o,x 1 ; 



n=2 



Figure 1: The plaquette presentation of the twisted boundary conditions. 



where the angular variables are changed in the following ranges 0, \ € [— 7r; +7r), € [0; 7r/2). 

The Abelian part of the lattice variables is represented by the diagonal components of the matrix 
and the condensate Abelian magnetic field influences the field 6> M (x), only. 

The second important task is to incorporate the magnetic flux in this formalism. The most natural 
way was proposed by 't Hooft [H]. In his approach, the constant homogeneous external flux <p in the 
third spatial direction can be introduced by applying the following t.b.c: 



It could be seen, the edge links in all directions are identified as usual periodic boundary conditions 
except for the links in the second spatial direction, for which the additional phase <p is added (Fig. 1). 
In the continuum limit, such t.b.c. settle the magnetic field with the potential A fl (x) — (0, 0, Hx 1 , 0). 
The magnetic flux <p is measured in angular units and can take a value from to 2ir. 

The lattice variables (in the Abelian parametrization) in the presence of the magnetic flux tp are 



where <Pfj,(x) — <p for the edge links at x = (xo,x%,L s , X3) with fj, = 2 and f^(x) = for other links. 

The total flux through the plane spanned by the plaquettes p, which affects the edge links at 
x = (xq, x%, L s ,xz) with // = 2, is 



U fl (L t ,x 1 ,X2,x 3 ) = U t _ l {0,x 1 ,x 2 ,x 3 ), 
U f _ l (x Q ,L s ,x 2 ,x 3 ) = U f _ l (x ,0,x 2 ,x 3 ), 



(8) 



Ufj,{xQ,xi,L s ,x 5 ) = e' i¥, [/ A1 (x ,a;i,0,a;3), 

U fl (xo,Xl,X2,L s ) = U fJ ,(xo,Xl,X2,0). 





(10) 



P = 0p U (x) = 6fj,(x) + Q v {x + ajl) - 9fj,(x + av) - 9 v {x). 



(11) 
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Eq. (|T0|) is the lattice analogue of the flux in the continuum: 

$ c = / ePtr^Fp,. (12) 
Js 

In this approach the variable tp describes a flux through the whole lattice plane, not just through an 
elementary plaquette. 

The t.b.c. for the components |(9|), 

Ufa) = cos((9 M (x) + tp ti (x))cos<pfa), (13) 
Ufa) = sin (j> fa) cob X fa), 
Ufa) = s'm(f>fa)smxfa), 



Ufa) = sin(0fa) + ipfa))coa<t>fa), 



read 



f Ufa)cosip - Ufa) simp for x = {xo,xi,L s ,x 3 } and fi = 2, 

U^iix) = \ (14) 
[ U°(x) for other links, 

„ f C/"(.T)sin^ + L73(a:)cos^ for x = {x , xi, L s , x 3 } and = 2, 

u£(x) = < (15) 
[ U^(x) for other links. 

The relations (fl4|) and have been implemented into the kernel of the MC procedure in order 
to produce the configurations with the magnetic flux tp. In this case the flux tp is accounted for in 
obtaining a Boltzmann ensemble at each MC iteration. 



3 Description of simulations and data fits 

The MC simulations are carried out by means of the heat bath method. The lattices 2 x 8 3 , 2 x 16 3 
and 4 x 8 3 at j3 = 3.0, 5.0 are considered. These values of the coupling constant correspond to 
the deconfinement phase and perturbative regime. To thermalize the system, 200-500 iterations are 
fulfilled. At each working iteration, the plaquette value ((5]) is averaged over the whole lattice leading to 
the Wilson action Q. Then the effective action is calculated by averaging over the 1000-5000 working 
iterations. By setting a set of magnetic fluxes tp in the MC simulations we obtain the corresponding 
set of values of the effective action. The value of the condensed magnetic flux ip m i n is obtained as the 
result of the minimization of the free energy density ((3J) • 

The spontaneous generation of magnetic field is the effect of order ~ g 0J. The results of MC 
simulations show the comparably large dispersion. So, the large amount of the MC data is collected 
and the standard x 2 -method for the analysis of data is applied to determine the effect. We consider 
the results of the MC simulations as observed "experimental data". 

The effective action depends smoothly on the flux tp in the region tp ~ 0. Therefore, the free energy 
density can be fitted by the quadratic function of the flux tp, 

F(tp) = F min +b(tp - tpmin) 2 - (16) 
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This choice is motivated also by the results obtained already in continuum field theory [13] where 
it was determined that free energy has a global minimum at tp ^ 0. The parametrization (fl6| is the 
most reasonable in this case. It is based on the effective action accounting for the one-loop plus daisy 
diagrams [13], 



having g 2 and (g 2 ) 3 ' 2 orders in coupling constant. Here, H is field strength (flux tp ~ H), T is the 
temperature, fi is the normalization point, IIoo(O) is the zero-zero component of the gluon polarazation 
operator calculated in the external field at the finite temperature and taken at zero momentum. The 
value of — 3, which was used, corresponds to a deep perturbation regime. So, a comparison with 
perturbation results is reasonable. The systematic errors in fitting function lfT6|) could come from not 
taking into account the high-order diagrams in lfT7j) . However, as it is well known [15J, the lack of an 
expansion parameter at finite temperature starts from the three-loop diagram contributions that is 
of g 6 order and could not remove an effect derived in g 2 and g 3 orders. As the finite-size effects are 
concerned, in the present investigation we just made calculations for two lattices 2 x 8 3 and 2 x 16 3 and 
have derived the same results for the tp m in (as it will be seen below). A more detailed investigation 
of this issue requires much more computer resources, which were limited. 

There are 3 unknown parameters, F m i n , b and <p m in in Eq. (|16|) . The parameter p m in denotes 
the minimum position of the free energy, whereas the F m i n and b are the free energy density at the 
minimum and the curvature of the free energy function, correspondingly. 

The value <p m in is obtained as the result of the minimization of the x 2 -function 

■Zfl? h \ {F m in + b(<fi — fmin) 2 ~ F(lfi)) 2 
X K- t min,0,<Pmin) = / , D(F(p-)) ' ^ ' 

where tpi is the array of the set fluxes and D(F(tpi)) is the data dispersion. It can be obtained by 
collecting the data into the bins (as a function of flux) , 

D(F( W )) = £ ^- P f\ (19) 

. ,_ , . Ubin ± 

where nu n is the number of points in the considered bin, Fu n is the mean value of free energy density 
in the considered bin. As it is determined in the data analysis, the dispersion is independent of the 
magnetic flux values ip. The deviation of (p m in from zero indicates the presence of spontaneously 
generated field. 

The fit results are given in the Table 1. As one can see, (p m in demonstrates the 2<r-deviation from 
zero. The dependence of (p m in on the temperature is also in accordance with the results known in 
perturbation theory: the increase in temperature results in the increase of the field strength [4]. 

The fits for the lattices 2 x 8 3 and 2 x 16 3 at j3 = 3.0 are shown in Fig. 2. The maximum- likelihood 
estimate of F(ip) by the whole data set is shown as the solid curve. In addition, all ip values are 
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Figure 2: x 2 -fit of the free energy density on lattice 2 x 8 3 (left figure) and 2 x 16 3 (right figure) 
for = 3.0 (grey regions describe the cp mi „ = 0.019+0 °" and cp mi „ = 0.0069^;^, at the 95% C.L., 
correspondingly) . 

divided into 15 bins. The mean values and the 95% confidence intervals are presented as points for 
each bin. The first 9 bins contain about 600-2000 points per bin. The large number of points in the 
bins allow to find the free energy F with the accuracy which substantively exceeds the dispersion, 
y/ D(F(ipi)) ~ 10 . It makes possible to detect the effect of interest. As it is also seen, the maximum- 
likelihood estimate of F(<p) is in a good accordance with the bins pointed, because the solid line is 
located in the 95% confidence intervals of all bins. 

The 95% C.L. area of the parameters F m i n (b for the right figure) and fmin is represented in 
Fig. 3. The black cross marks the position of the maximum-likelihood values of F m , n (b for the right 
figure) It can be seen that the flux is positive determined. The 95% C.L. area becomes 

more symmetric with the center at the F m i n , b and f m i n when the statistics is increasing. This also 
confirms the results of the fitting. 



4 Discussion 

The main conclusion from the results obtained is that the spontaneously created temperature-depen- 
dent chromomagnetic field is present in the deconfinement phase of QCD. This supports the results 
derived already in the continuum quantum field theory [H [12] and in lattice data analysis |9j . 

Let us first discuss the stability of the magnetic field at high temperature. It was observed in Refs. 



Table 1: The values of the generated fluxes (pmm for different lattices (at the 95% C.L.). 





2 x 8 3 


2 x 16 3 


4 x 8 3 


= 3.0 


019+ 0013 
u.u±»_ 012 


69+ 0022 
u.uuoy_ 0057 


no5 +0 - 005 

U.UUO_ 003 


= 5.0 


020+ 011 
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Figure 3: The 95% C.L. area for the parameters F min and <p m i n , determining the free energy density 
dependence on the flux (p m in on lattice 2 x 8 3 for j3 = 3.0 (left figure). The 95% C.L. area for the 
parameters F m j n and b , determining the free energy density dependence on the flux <-p m in on lattice 
2 x 16 3 for (3 = 3.0 (right figure). 

[H [12] that the stabilization happens due to the gluon magnetic mass calculated from the one-loop 
polarization operator in the field at temperature. This mass has the order rn 2 nagn ~ g 2 (gH) 1 / 2 T ~ 
g 4 T 2 as it should be because the chromomagnetic field is of order (gH) 1 / 2 ~ g 2 T [4]. The stabilization 
is a nontrivial fact that, in principle, could be changed when the higher order Feynman diagrams to 
be accounted for. Now we see that the stabilization of the field really takes place. 

Our approach based on the joining of calculation of the free energy functional and the consequent 
statistical analysis of its minimum positions at various temperatures and flux values. This overcomes 
the difficulties peculiar to the description of the field on a lattice. Here we mean that the field strength 
on a lattice is quantized and therefore a nontrivial tuning of the coupling constant, temperature and 
field strength values has to be done in order to determine the spontaneously created magnetic field. 

We also would like to note that in the present paper the flux dependence on temperature remains 
not investigated in details. This is because of the small lattice size considered. That restricts the 
number of points permissible to study. However, at this stage we have determined the effect of 
interest as a whole. Even at the small lattice, one needs to take into consideration thousands points of 
free energy (that corresponds to an analysis of 5-10 millions MC configurations for different lattices) to 
determine the flux value >p m in at the 95% C.L. In case of larger lattices this number and corresponding 
computer resources should be increased considerably. This problem is left for the future. 

As we mentioned in Section 2, the finite-size effects have not been investigated in detail. However, 
these effects are important near the phase transition temperature. They make difficult to distinguish 
a first-order phase transition from a second-order one. In our case, the temperature is far from T c . 
The fact that external field penetrates the Coulomb phase is well known |10[ [T6] , so the only really 
new thing is that this field is spontaneously created. It was first observed in continuum [4], where 
the field strength of order g A in coupling constant was established. Finite-size effects are not able to 
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remove this result. The values (p m in obtained on the lattices 2 x 8 3 and 2 x 16 3 (see the Table 1) are 
in a good agrement with each other, within the statistical errors at 95% C.L. 

One could speculate that the lattice sizes 2 x 8 3 and 2 x 16 3 are not sufficient. However, these 
lattice sizes were used in the Refs. [17]. The main aim of present paper is to show a possibility of 
spontaneous generation of chromomagnetic field at high temperature in lattice simulation, which was 
investigated already by perturbative methods [4j [8]. 

It is interesting to compare our results with that of in Ref. [10] where the response of the vacuum 
on the external field was investigated. These authors have observed in lattice simulations for the 
SU(2)- and S'C/(3)-gluodynamics that the external field is completely screened by the vacuum at low 
temperatures, as it should be in the confinement phase. With the temperature increase, the field 
penetrates into the vacuum and, moreover, increase in temperature results in existing more strong 
external fields in the vacuum. On the other hand, increase in the applied external field strength leads 
to the decreasing of the deconfinement temperature. These interesting properties are closely related 
to the studies in the present work. Actually, we have also investigated the vacuum properties as an 
external field problem when the field is described in terms of fluxes. This was the first step of the 
calculations. The next step was the statistical analysis of the minimum position of free energy, in 
order to determine the spontaneous creation of the field. In fact, at the first step we reproduced the 
results of Refs. [10] (in terms of fluxes). 

Note that the present investigations also correspond to the case of the early universe. They support 
our previous results on the magnetic field generation in the standard model [13J and in the minimal 
supersymmetric standard model [14]. As it was discussed by Pollock [5], the field generated by this 
mechanism at the Planck era might serve as a seed field to produce the present day magnetic fields 
in galaxies. 

We would like to conclude with the note that the deconfinement phase of gauge theories is a very 
interesting object to study. The temperature dependent magnetic fields, which are present in this 
state, influence various processes that should be taken into consideration to have an adequate concept 
about them. 
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